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ABSTRACT 

We present new high-resolution (~(y.'09) ii-band imaging observations of the circumstellar 
disk around the T Tauri star SU Aur. Our observations with Subaru-HiCIAO have revealed the 
presence of scattered light as close as 0'.'15 (~20 AU) to the star. Within our image, we identify 
bright emission associated with a disk with a minimum radius of ~90 AU, an inclination of ~35° 
from the plane of the sky, and an approximate P.A. of 15° for the major axis. We find a brightness 
asymmetry between the northern and southern sides of the disk due to a non-axisymmetric disk 
structure. We also identify a pair of asymmetric tail structures extending east and west from the 
disk. The western tail extends at least 2"5 (350 AU) from the star, and is probably associated 
with a reflection nebula previously observed at optical and near-IR wavelengths. The eastern 
tail extends at least 1"(140 AU) at the present signal-to-noise. These tails are likely due to an 
encounter with an unseen brown dwarf, but our results do not exclude the explanation that these 
tails are outflow cavities or jets. 

Subject headings: protoplanetary disks — stars: individual (SU Aur) — stars: pre-main sequence — 
polarization 
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1. INTRODUCTION 

Direct images of scattered light from protoplan¬ 
etary disks provide valuable information about 
their surface geometries. Advances in instru¬ 
mentation have produced high angular resolution, 
high sensitivity observations, revealing a variety of 
substructures in protoplanet ary disks associated 
with pre-main sequence (PMS) stars. In particu¬ 
lar, near-infrared (IR) imaging polarimetry with 
adaptive optics has been used extensively dur¬ 
ing the Strategic Explorations of Exoplanets and 
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Disks with Subaru (SEEDS) survey with HiCIAO 


(Tamura 2009). This survey has been successful in 


observing diverse disk morphologies, from uniform 
(e.g., MWC 480— Kusakabe et al.||2012 ) to com¬ 
plex, including spirals (e.g., AB Aur, SAO 206262, 
MWC 758— Hashimoto et al. 2011 |Muto et al. 


20l2l|Grady et al.||2013[ ), gaps or rings (e.g., PDS 


70, 2MASS J1604, Oph IRS 48— Hashimoto et al.| 


2012 Mayamaet al. 2012 Follette et al.|2015 ), and 
other asymmetric flux distributions (e.g., UX Tau 
A, SR 21, RY Tau, Oph IRS 48— [Tanii et al.|2012 


Follette et al. 2013[ Takami et al. 1120*13 Follette 


et al.||2015 ). Some of these structures may be due 
to dynamical disk-planet interactions providing 


hints of ongoing planet formation (e.g., Hashimoto 

et al. 2011 

2012 

Muto et al. 2012J Mayama et al. 

2012 

Grady et al. 2013 

)• 


SU Aur is a classical T Tauri star located in 
the Taurus-Auriga star formation region at a dis¬ 
tance of ~140 pc (Bertout and Genova 20061. This 


star-disk system shows large IR and UV excesses 
indicating the existence of a circumstellar disk and 
ongoing active mass accretion, respectively (e.g., 


Bertout et al.|1988;|Gaivet et al. 2004; Jeffers et al. 

2014). Table 

properties. 

SU Aur is 

|]J summarizes its stellar and disk 

also associated with a bright re- 

flection nebula (Nakajima and Golimowski 

1995 

Grady et al.| 

2001 1 Chakraborty and Ge 

2004 

Jeffers et al. 

20141. Nakajima and Golimowski 


(1995) revealed an associated cometary nebulos¬ 


ity at optical wavelengths, extending to the west 
and southwest at small (<10") and large distances 
(10"-30"), respectively. Optical high-resolution 
observations using coronagraphy or imaging po- 


larimetry 

(Grady et al. 

2001 

Chakraborty and 

Ge 2004 

offers et al. 201 

4) have resolved the ex- 


tended nebulosity toward the west within 5" of 
the star. Possible origins of the observed struc- 


Chakraborty and Ge 2004), part of the circmustel- 

lar disk (Chakraborty and Ge 2004 

, or part of a 

remnant envelope or parent cloud 

(Jeffers et al. 

2014). Chakraborty and Ge ( 

2004 

also applied 


point-spread function (PSF) subtraction to seeing- 
limited near-IR images, and report the presence of 
an extended nebulosity to the southwest, close to 
the star. 

In this paper, we present near-IR imaging po¬ 
larimetry of SU Aur. Using Subaru-HiCIAO, we 
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have revealed scattered light from the disk, as well 
as two tails which may be due to tidal interaction 
with an unseen brown dwarf. The rest of the paper 
is organized as follows. In section 2, we summarize 
our observations and data reduction. In section 
3, we show the observed Polarized Intensity (PI) 
flux distribution and polarization. We highlight 
the discovery of salient features including the disk 
and tail structures. In section 4, we discuss the 
orientation of the disk and the nature of the tail 
structure. 


2. OBSERVATIONS AND DATA RE¬ 
DUCTION 

SU Aur was observed in //-band (1.65 /iin) 
on 2014 January 19 UT using the high-contrast 
imager HiCIAO. The adaptive optics system 
(A0188) provided a stable stellar PSF (FWHM = 
0708, Strehl ratio = 0.2-0.3). We used a combina¬ 
tion of polarimetric and angular differential imag¬ 
ing (PDI+ADI) modes. A dual Wollaston prism 
was used to split incident light into two image 
pairs that are orthogonally and linearly polarized, 
each with a 5" x5" field of view and a pixel scale 
of 9.5 mas pixel -1 . As for several other SEEDS 
observations, linear polarization was measured by 
rotating the half-wave plate to four angular posi¬ 
tions: 0°, 22.5°, 45°, and 67.5°. We obtained 12 
full waveplate cycles, taking a 30-sec exposure per 
waveplate position, using a 10% neutral density 
(ND) filter. No coronagraphic mask was used in 
order to image the innermost region around the 
central star. Four additional sets of waveplate 
cycles were obtained using a 1% ND filter to mea¬ 
sure the PSF FWHM and the total Stokes / flux 
density of the star (/*). 

The data were reduced using the standard 


method for ADI+PDI (Hinkley et al. 20091. 


The measured Stokes /, Q, and U parameters 
were used to compute the more accurate intrin¬ 
sic Stokes parameters. We then calculated the 
polarized intensity (PI = yj ( Q 2 + U 2 )), polariza¬ 
tion (PI//) and the angle of polarization (dp = 
0.5xarctan([// Q )). The instrumental polarization 
of HiCIAO at the Nasmyth platform was corrected 


following Joos et al. (2008) with errors of <0.1%. 


To increase the signal-to-noise, we convolved 
the Q and U images with a Gaussian having a 
FWHM of 0'.'04. As a result, the angular resolu¬ 


tion of the images shown in later sections is 0 / /09, 
slightly larger than the instrumental resolution. 


We also applied the LOCI algorithm ( Lafreniere| 


et al. 2007) to the Stokes / images to search for 


faint companions, with parameters as in Table [2] 
The image rotation during the observations was 
small (~8°) hence we needed a small N a (the min¬ 
imum displacement distance in the azimuthal di¬ 
rection) to search for a companion close to the 
star (down to r=0'.'17 and 0736 for V CT =0.1 and 
0.5 x FWHM of the observations, respectively). 
The other LOCI parameters were the same for 
the four sets and were standard and optimum 
for SEEDS observations. We embedded artificial 
point sources at different radii and derived upper 
limits for 3-u fluxes of the companion of mfl-=~15, 
~19, 21.0 and 21.5 at r=0"16, 075, 1", and 2" from 
the star, respectively. The detection limits were 
almost identical between different parameter sets. 


3. RESULTS 


Figure |T] shows the observed PI flux distribu¬ 
tion overlaid with polarization vectors. There is 
bright emission at the center with a radius of ~100 
AU. This emission probably originated from a 
protoplanetary disk elongated to the north-south, 
slightly offset in the counter-clockwise direction. 
The peak emission in the disk is located at r~0720 
(^28 AU) north of the star and has PI//* of 
2.9xl0 -6 pix -1 . This translates to surface bright¬ 
ness in PI of ~80 mJy arcsec -2 . This is compara¬ 
ble to the brightest disks in the SEEDS observa¬ 


tions summarized in Takami et al. (20141. More¬ 


over, the distribution shows a deficit of emission 
perpendicular to the semi-major axis similar to PI 
distributions observed in other disks with inter¬ 


mediate inclination angles (e.g., 

Hashimoto et al. 

2012 

Kusakabe et al. 2012 

Tanii et al. 2012). This 


deficit likely arises because along the minor axis, 
the light is dominated by forward- and backward- 
scattered photons, which are intrinsically less po¬ 
larized (e.g., Takami et al.|2013 ), and does not re¬ 
quire an actual deficit of material in the disk. We 
detect a brightness asymmetry between the disk’s 
northern and southern sides: the northern side is 
brighter than the southern side by a factor of 2. 


To analyze the inclination and orientation of 
SU Aur’s disk, we fit an ellipse to the outermost 
contour (PI//* = 10 -7 pix -1 ) by a simple visual 
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analysis. This fit yields angular separations of 177 
and 174 (~240 AU and ~200 AU) for the major 
and minor axes, respectively, with a position an¬ 
gle (P-A.) of ~15° for the major axis. The ratio 
between the minor and major axes indicates an ap¬ 
proximate disk inclination of 35°. We do not make 
a more quantitative or accurate measurement for 
the disk orientation or inclination because of the 
asymmetric geometry seen in the contour. 

We extracted the radial PI//* profiles wherein 
we averaged a swath of 11 pixels (~0 / /l) along the 
disk’s major axis on both sides of the star. A 
power-law function was fitted to each profile, as 
shown in Figure [2] Between 0725 and 076, we 
find the following power-law indices that best fit 
the observed PI profiles: ~-3 and ~-2 at P.A. of 
15° and 195°, respectively. The indices we mea¬ 
sured are roughly consistent with SEEDS studies 
of several disks (SAO 206462, MWC 480, MWC 


758 — Muto et al.||2012) |Kusakabe et al.||2012| 


Grady et al. 2013) as well as some other obser¬ 


vations of disks associated with Herbig Ae and Be 
stars (Fukagawa et al. 2010| . Compared with P.A. 
= 15°, the observed profile shows a larger devia¬ 
tion from the power-law fit at P.A.=195°(i.e., shal¬ 
lower and steeper slopes at inner and outer radii, 
respectively). A different brightness distribution 
between these two P.A.s indicates that the disk 
structure is non-axisymmetric. 


Figure [2] also shows a drop-off in the observed 
flux at 076-077 from the star. This indicates 
a lower limit for the disk radius of ^90 AU at 
a distance of 140 pc. This is consistent with 
disk models for millimeter emission and optical-to- 
millimeter spectral energy distributions by |Ricci| 
etak1 ( [M0| (100-300 AU) and | Jeffers etak] poTIj ) 


(500 AU), respectively. Note that the absence of 
flux in the outer region does not imply the absence 
of disk material as it may be due to self-shadowing 
of the disk ( |Takami et al.|[2014| . 

In Figure]!] a long, faint tail (PI//* < 0.5x 10~ 7 
or 1.4 mJy arcsec -2 ) extending westward at least 
275 (350 AU) can also be seen. We also find an 
extended component that is marginally detected 
from the north-east to the east side of the disk. 
The tail in the west is likely associated with the 
extended reflection nebula previously observed at 


optical and near-IR wavelengths I 

Nakajima and 

Golimowski|1995j Grady et al. 2001[ Chakraborty 

and Ge 2004 Jeffers et al. 2014). 

However, the 


distribution of emission appears significantly nar¬ 
rower in the north-south direction than previous 
observations in the optical and near-IR. 


To better display structures in the faint outer 
regions, we convolve the PI image with a Gaussian 
of FWHM=07l and then scale the PI flux at each 
pixel by R 2 , where R is the projected distance 
from the star (Figure [3]). The figure clearly shows 
two tails extending east and west. In Figure[3j the 
tails appear to be associated with the edge of the 
disk emission, however, observations with higher 
signal-to-noise are required to confirm this. 


Figure [l] also shows the polarization vector 
map. The vectors in the disk show a centrosym- 
metric pattern, as has been observed in several 
other disks with near-IR imaging polarimetry 
(e.g., Hashimoto et al.||2011 |. This indicates that 
the disk is seen in scattered light that originated 
from the central star. Similarly, the vectors in the 
tail indicate that it is illuminated by the star-disk 
system. The disk has lower fractional polariza- 
tion (<10%) than the western tail (up to ^30%). 
This is because the Stokes / flux within 1" of the 
star is dominated by the star’s PSF halo, lowering 
the degree of polarization (PI//). Hence, the vec¬ 
tors represent the lower limit of the true degree of 
polarization. 


4. DISCUSSION 


Jeffers et al. (20141 conducted optical imaging 


polarimetry of SU Aur and found a disk-like mor¬ 
phology elongated in the east-west direction at a 
larger angular scale (radius of ~2"). The disk-like 
structure within l"of the star seen in the near-IR 
(Figure [l]) is almost perpendicular to this struc¬ 
ture in the optical. As the dust opacity is higher 
in the optical than in the near-IR, the disk may 
be partially obscured by dust in the envelope at 
optical wavelengths. Thus, near-IR observations 
should have an advantage in observing scattered 


light close to the disk surface (Fukagawa et al. 


2004). This explanation is supported by the fact 


that the power-law indices of the radial PI dis¬ 
tribution in the near-IR are similar to other disk 
systems (Section 3). 


Tails like those in Figures |T] and [3] are not usu¬ 
ally observed in scattered light associated with 
young stellar objects (YSOs). An exception may 
be a jet-like feature in the Z CMa system observed 
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by Millan-Gabet and Monnier (2002), which the 


authors attributed to a cavity wall. In the re¬ 
maining section we discuss possible origins of these 
tails: (1) outflow cavities, (2) collimated jets, and 
(3) tidal tails due to a (sub-)stellar encounter. 

Before discussing the nature of the SU Aur tails 
in detail, we estimate the dust mass of the western 


tail. Takami et al. (20131 derived the following 


equation to estimate the dust mass observed in 
scattered light using imaging polarimetry: 


f npi{ r) r 2 / PI\~ 1 

TTldust = / ( ) 

J 71* /^ext ' 1-0 ' 


dr 


Cl) 


where rridust is the dust mass; npi(r) is the num¬ 
ber of PI photons at each position; n* is the num¬ 
ber of stellar / photons; r is the distance to the 
star; n ext is the dust opacity; PI/Iq is the frac¬ 
tion of the PI flux normalized to the incident flux 
on the dust grains. Here we assume the interstel¬ 


lar du st size distribution measured by Kim et a l. 
(1994). This yields K ext = 5.3 x 10 3 cm 2 g _1 and 


PI/I 0 ~ 0.01 at a wavelength of 1.65 pm (Takami 


et al. 2013). Integrating Equation (1) from disk 
edge to the tail end in Figure [l] (~110 - 350 AU) 
yields a dust mass of ~6 xlO -8 M 0 , if we use 
the projected distance for r at the distance of 140 
pc. This mass is smaller than the disk’s dust mass 
inferred from millimeter interferometry (see Table 
[l]) by a factor of 2000-5000. However, like proto¬ 
planetary disks, a significantly larger mass might 
exist behind the scattering layer. The above dust 
mass should therefore be a lower limit for the tail. 

SU Aur is known to be associated with a re¬ 
flection nebula in the east-west direction, approx¬ 
imately the same direction as the tails in Figures 


lar scales ( 

Nakajima and Golimowski 

1995 Grady 

et al. 2001 

Clrakraborty and Ge[2004 

Jeffers et al. 

2014). Such nebulosity associated with YSOs is of- 


ten attributed to an outflow cavity (e.g., 

Tamura 

et al. 1991 

Hodapp 

1994 Lucas and Roc 

re 1998] 

Padgett et al. 

1999 

). However, this explanation 


may not apply to these tails for several reasons. 
First, the tails shown in Figures [T| and [3] are sig¬ 
nificantly narrower than other YSO outflow cav¬ 
ities, which generally have a wide opening angle. 
Secondly, an outflow cavity in scattered light is 
usually observed toward a single direction unless 


we see the disk edge-on (e.g., Fischer et al. 1994 


Whitney et al. 2003). However, the tails in Fig¬ 


ures [T] and [3] are seen in two opposite directions 
despite the fact that the disk has an intermediate 
viewing angle. 


Another possible explanation for the tails may 
be that these are collimated jets. Many YSOs 
are associated with such jets (e.g., Frank et al. 
2014 for a review), however, these are not usu¬ 


ally observed in scattered light. To further inves¬ 
tigate this explanation, we estimate the mass ejec¬ 
tion rate and compare it with the disk accretion 
rate. Using Equation (1) we derive a dust mass of 
3xl0 23 g AU -1 along the western tail at 200 AU 
from the star. This would yield a mass ejection 
rate of 2xl0~' yr _1 , assuming a typical jet 

(e.g., Hartigan et al. 1995 


velocity of 50 km s 


-1 


Hirth et al. 1997) and gas-to-dust mass ratio of 


100. This is 7-40 times larger than the disk ac¬ 
cretion rate measured by Calvet et al. (2004) and 
Ricci et al. (2010) (see Table|lJ). The mass ejection 


to mass accretion ratio would be even larger if the 
tail contains more gas and dust than is inferred 
from the scattered light. However, such a high 
mass ejection to mass accretion ratio is not feasi¬ 
ble for PMS stars, which usually show a mass ejec¬ 
tion to mass accretion ratio of ^0.1 (e.g., Calvet 


1997). Therefore, jets are not likely to explain the 


tail structure either, unless the jet velocity is sig¬ 
nificantly lower (e.g., <10 km s _1 ) than for other 
similar objects. 


Tidal interaction with a star or a brown dwarf 
may alternatively explain the observed tail struc¬ 
tures. Kinematic simulations of stellar encounters 
show that this physical process produces tails sim¬ 
ilar to those in Figures[l]and[3](c.g, Pfalzner 2003 


Forgan and Rice 2009). A tidal tail was observed 


in the RW Aur disk system using millimeter inter¬ 
ferometry (Cabrit et al. 2006). A difficulty with 


this interpretation is that there is yet no observa¬ 
tional evidence of another star within 30" of SU 


Aur (Nakajima and GolimowskiJ 1995 Grady et al. 


|200l( Clrakraborty and Ge|2004j| Jeffers et al.[2014j 

see also data archive for 2MASS). However, an en¬ 
counter with an unseen brown dwarf (BD) could 
produce tails observable in scattered light. We 
estimate an upper limit BD mass of 0.6-2xl0” 3 
Mq for our field of view of observations with a de¬ 
tection limit of 21.0-21.5 mag (Section 2) and the 
COND models for 1-10 Myr (Baraffe et al. (2003). 
The same models and the completeness limit of the 
2MASS Point Source Catalog (15.1 mag) yield an 
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upper limit mass of 0.8-2x10 2 M @ . 

A (sub-)stellar encounter would affect the disk 
structure, and therefore affect mass accretion and 
planet formation (e.g, Pfalzner 2003 Forgan and 
Rice] [20091 [SOToj |Breslau et al.||2014[ ). SU Aur 
would be a good, rare observational example to 
test these theories in detail. Confirmation of this 
physical process in this system would require deep 
imaging to search for a nearby brown dwarf, and 
kinematic information for the disk and the tail 
from millimeter interferometry. A search for jet 
emission at optical and near-IR wavelengths would 


also be a good test for the jet hypothesis (Frank 
et~aJ7|[2014l ) . 


5. CONCLUSIONS 


We present H-band imaging polarimetry of 
scattered light from the circumstellar disk around 
the T Tauri star SU Aur. Our high-resolution 
(0'.'09) image of polarized intensity has allowed us 
to observe scattered light from the disk as close 
as 0715 (~20 AU) to the star. We identify bright 
emission probably associated with the disk, elon¬ 
gated along an approximate P.A. of 15°. From its 
distribution, we estimate a minimum disk radius 
of ~90 AU and a disk inclination of ^35°. As 
observed in several protoplanetary disks, the sur¬ 
face brightness along the disk’s major axis is ap¬ 
proximately proportional to r~ 2 to r -3 at r= 072- 
0'/65, where r is the projected distance to the star. 
The brightness distribution of the scattered light 
is asymmetric between the disk’s northern and 
southern sides indicating a non-axisymmetric disk 
structure. 

Our near-IR images also show two tails extend¬ 
ing to the east and west of the disk. The western 
tail extends at least 2 , . , 5 (350 AU) from the star, 
and is probably associated with an extended re¬ 
flection nebula previously observed at optical and 
near-IR wavelengths. The eastern tail extends at 
least 1"(140 AU) at the present signal-to-noise. 
We estimate a lower limit of the dust mass of the 
western tail (r=110-350 AU) of ~6 x 10~ 8 Mq 
using an interstellar dust model and assuming 
that the tail lies close to the plane of the sky. 
This mass is significantly lower than the disk’s 
dust mass inferred from millimeter interferometry 
(1 — 3 x 10~ 4 Mq). 

Possible explanations for the tails include (1) 


outflow cavities, (2) collimated jets, and (3) tidal 
interaction with a brown dwarf. The morphol¬ 
ogy and brightness favor the last explanation, but 
these do not rule out the other possibilities. Con¬ 
firmation requires searching for a nearby brown 
dwarf with deep imaging, and observing the kine¬ 
matics of the disk and tails using millimeter inter¬ 
ferometry. 
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Fig. 1.— PI image of SU Aur in log-scale overlaid with polarization vector map. The PI flux at each pixel is 
scaled by the spatially integrated stellar I flux (PI//*=10 - ' corresponds to 2.8 mJy arcsec -2 ). The contours 
are shown for a log scale of PI//* of -7.0, -6.75, -6.5 and -6.25. The dotted white ellipse is adjusted to 
the outermost contour (see text for details). The region within r =0715 (~20 AU) is masked to show only 
regions where the measurements are reliable. 
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Fig. 2.— Radial PI//* profiles along the semi¬ 
major axis of the disk. The filled circles and solid 
squares show the measured value at P.A.s of 15° 
and 195°, respectively. The solid and dashed lines 
show power-law fits made at r=0 , . , 25-0 , .'6, where r 
is the projected distance to the star. 



-1.0 -0.5 0.0 0.5 1.0 


R.A. Offset (arcsec) 

Fig. 3.— The outer disk and tails, highlighted 
by scaling the PI flux at each pixel by /? 2 , where 
R is the projected distance to the star. The im¬ 
age is shown within 1"3 of the star, where this 
image processing does not significantly enhance 
background fluctuation. The color scale is arbi¬ 
trary. The white contours show the original PI 
flux distribution shown in Figure [T] 
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Table 1: Parameters for SU Aur 


Parameter 


Reference 

Spectral Type 

G2III 

1 

Stellar Mass 

1.9 ±0.1 M q 

2 

Age 

6.8 ±0.1 Myr 

2 

Dust mass of the disk 

1-3 xlO" 4 Mq 

3 

Mass accretion rate 

0.5-3xl0“ 8 Mq yr" 1 

4,5 

Distance 

143^3 pc 

6 

Ay 

0.9 mag 

2 

H mag. 

6.56 mag 

7 


Referen ces. — [1] |HerbIg| (|1952|) ; | 2] Bert out et al.| ( |2007[ ) ; [3] |Ricci et al.| < |2010[ ) ; [4] |Calvet et al.| ( |2Q04[ ) ; [5] | Ricci et al.| 
( |2010| ); [6] |Bertout and Genovaj ( |20Q6p ; [7] 2MASS All Sky Catalog of Point Sources 


Table 2: LOCI Parameters 


Parameter 3, 

Values b 

N a 

150 x FWHM 

dr 

5 x FWHM (r < 50 x FWHM) 


100x FWHM (r > 50 x FWHM) 

9 

1 

N a 

0.1/0.3/0.5/1.0 x FWHM 


a N A ••• the size of the optimization area; dr ... the width of the subtraction area; g ... the ratio of radial and azimuthal width 
for the optimization area; N& ... the minimum displacement distance in the azimuthal direction ( [Lafreniere et al.|2007| >. 
b FWMH is that of the observations (O'/OS). 
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